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oscillator with respect to the classical Reissner-Nordstrom black hole. 
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I. INTRODUCTION 



Quantum theory and General Relativity are two cornerstones of modern physics that for 
more than a century have contributed to increase our knowledge of the Universe as never 
before in the human history With the help of the Quantum Theory we can explain micro- 
world phenomena, and the General Theory of Relativity allows us a deep understanding 
of the Universe at cosmological scales. Unfortunately, this two beautiful theories resist 
all attempts to bring them together. A unified theory of gravity and the quantum world 
would be very important to describe, for example, the origin of the Universe and its later 
development. 

There are other simple phenomena that can be very interesting to describe by the future 
Quantum Gravity. Among other things, from the classical side, it is well known that the 
response of a black hole to small perturbations at intermediate times is characterized, under 
suitable boundary conditions, by a discrete set of complex frequencies called quasinormal 
frequencies, that depend only upon the parameters of the black hole [l -^J. From the quan- 
tum side, it would be interesting to see what changes appear in the evolution of quantum 
black holes under perturbations. Specially interesting is the behavior at intermediate times 
dominated by quasinormal response, because apart from allowing to us to gain some valu- 
able information about these objects, the quasinormal spectrum permits investigation of 
the black hole stability against small perturbations. Several numerical methods have been 
developed to study such interesting problem |5K7|]. 

Quasinormal modes appears to be important in other contexts, as for example, the 
AdS/CFT correspondence, where the inverse of imaginary part of quasinormal frequencies 
of AdS black holes can be interpreted as the dual CFT relaxation time js] {9]. 

In a previous paper we considered the influence of vacuum polarizations effects due to the 
backreaction of a quantum massive scalar field of large mass upon the quasinormal modes 
of elecrically charged black hole solutions obtained solving the semiclassical Einstein field 
equatio ns, with the quantum renormalized stress tensor of the quantized matter field as a 



source [10|. Such an influence appears essentially as an appreciable shift in the quasinormal 
frequencies that decreases as the bare black hole mass increases, and that not have a strong 
dependance upon the quantum field parameters, leading to the conclusion that the quantum 
corrected black holes are less oscillatory with respect to its classical counterparts. Another 



previous work along similar lines was done by Konoplya [ll|], for the BTZ black hole dressed 
by a massless scalar field, but in this case he considered the influence of particle creation 
around the event horizon, an effect that dominates over the vacuum polarization effect for 
massless fields. 

To solve the backreaction problem in semiclassical gravity, we need to know the functional 
dependence of the renormalized stress energy tensor of the quantum field surrounding the 
classical compact object on a wide class of metrics |l2|]. Unfortunately, this is a very difficult 
problem, and up to now, there exist only approximate methods to develop a tractable 
expression for this quantity 12H19I]. Since the pioneering work of York [20j, who solved 



the semiclassical Einstein equations for a Schwarzschild black hole dressed by a massless 
conformally coupled scalar field, using for the quantum stress energy tensor the results 
given earlier by Page 13J, there are some related works in the literature, both for massless 
and massive quantum fields of different values of the spin parameter. To see the effects 
of the backreaction upon the black hole response to small perturbations, quantum massless 
fields as sources of the quantum corrections are not the most suitable candidates, because the 
semiclassical metric components diverge as r — > oo and to obtain the correct solutions to the 
backreaction problem we need to impose some sort of boundary to the system under study, 
a feature that causes a change the quasinormal spectrum. A different situation happens in 
the case of very massive fields, for which the vacuum polarization effects are not difficult to 
compute constructing the quantum stress energy tensor by means of the Schwinger-DeWitt 
expansion of the quantum effective action, whenever the Compton's wavelenght of the field 
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is less than the characteristic radius of curvature 

It is important to mention that in semiclassical gravity the unavoidable effects due to 
the metric fluctuations and the associated graviton contributions to the complete quantum 
stress-energy tensor are ignored, a fact that is usually justified considering that there exists 
a regime in which the gravitational field can be regarded as a classical entity, and the effects 
of the remaining matter fields after quantization can be taken as quantum corrections to the 
bare metric. Using the quantum stress-energy tensor of the matter fields as a perturbation 
in the right hand side of the semiclassical Einstein equations, we can obtain a perturbative 
solution to the backreaction problem up to first order, and determine what changes appear 
in some important quantities as the mass, the location of the event horizon and the Hawking 
temperature of the quantum corrected solution. 



In this paper we study the effects that vacuum polarization of very massive scalar, vector 
and spinor fields cause on quasinormal modes of quantum corrected Reissner-Nordstrom 



black holes in four dimensions. This is the sequel of our previous work [10| in which we 
focus on the quantum scalar field case. In the first section we review the Schwinger-DeWitt 
technique to obtain the one-loop approximation for the effective action for massive fields 
in the large mass limit, and present the particular results obtained for a classical Reissner- 
Nordstrom black hole background. In section II we solve the backreaction problem to obtain 
the metric that describes the spacetime geometry of an electrically charged semiclassical 
black hole. Section III is devoted to the calculation of the massless test scalar quasinormal 
frequencies in this semiclassical background, by sixth order WKB method. Finally in Section 
IV we give the concluding remarks and comment on related problems to be studied. 

In the following we use for the Riemann tensor, its contractions, and the covariant deriva- 
tives the sign conventions of Misner, Thorne and Wheeler [24|. Our units are such that 

n = c = G = i. 

II. RENORMALIZED STRESS ENERGY TENSOR FOR QUANTUM MASSIVE 
FIELDS 

In the following we consider the quantization of massive scalar, vector and spinor fields 
in the l arg e mass limit. The results for the massive scalar field can be found in our previous 



works 
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15| . and for this reason we will be concerned only with the vector and spinor 
cases. The action for a single massive vector field with mass m v in some generic curved 
spacetime in four dimensions is 

Sv = -J d^V^g V + \ m l A ^ ■ (!) 

The equation of motion for the field have the form 

C(V)A M = , (2) 
where the second order operator V£ (V) is given by 

V»{y) = ^U-^ v W-W v -m 2 X , (3) 
where □ = g^V '„V 1 v is the covariant D'Alembert operator, V M is the covariant derivative. 
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For single massive neutral spinor field the action is: 

S f = % -J d A x^~g4> [7^0 + m f <j>] . (4) 

In the above expression, provides a spin representation of the vierbein group and = 0*7, 
where * means transpose. The Dirac matrices 7 and7^ satisfy the usual relation [7^, 7"], = 
2g^ u I, where I is the 4x4 unit matrix. 

The covariant derivative of any spinor ( obey the conmutation relations 23|, |25 | 
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v M v,c - v„v M c = gfffcfli* a V > ( 5 ) 
v„ v ff v„c - v CT v, v^c = ld M R a \^ u c + V P CR, p va , (6) 

+ V„V M C^' aT , (7) 
and so forth, where $[ a ,p] = 7 [7a>7 / a]_ are the generators of the vierbein group, [ , ]_ is 
the commutator bracket, and R a ° = h a a .h\R (rr with h a g the vierbein which satisfies 
h a ij,h a u = g^v The covariant derivatives of 7, 7^ and $\ a ,j3\ vanishes. The equation of motion 
for the field derived from the action reads 

(yV„ + m/)0 = O . (8) 

The operator Z)y that gives the evolution of the spinor function in flSj) is: 

£ / = 7 M V M + m / . (9) 

The usual formalism of Quantum Field Theory gives an expression for the effective action 
of the quantum fields Ap, as a perturbative expansion, 

r(^,0) = s(^,0) + ^r w (A /3 ,0) , (10) 

k>l 

where S(Ap,<p) is the classical action of the free fields. The one loop contribution of the 
fields A/3, to the effective action is expressed in terms of the operators (|3j) and (jHJ) as: 

r (1) = % - In (ptiV) + % - In (©etD) , (11) 

where QctF = exp(TrlnF) is the functional Berezin superdeterminant of the operator 
F, and TrF = (—lfF- = J d 4 x (— 1) A F£(x) is the functional supertrace [21]. If the 



Compton's wavelength of the field is less than the characteristic radius of spacetime curvature 
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26] , we can develope an expansion of the above effective action in powers of 
the inverse square mass of the field. This is known as the Schwinger-DeWitt approximation, 
and can be applied to "minimal" second order diferential operator of the general form 

i^(V) = ^D-m 2 <^ + Q£ , (12) 

where Q%(x) is some arbitrary matrix playing the role of the potential. 

Unfortunately, this is not the case of operators ([3]) and OH]). In the case of ([3]), the presence 
of the nondiagonal term turn it to be a nonmimal operator. 

By fortune we can put fl3]) as function of some minimal operators, if we note that it 
satisfies the identity VJf (V) (m 2 v 5£ - V„ V) = m\ - R£ - m 2 v 5^) . Then the one loop 
effective action for the nonminimal operator ([3]) omitting an inessential constant can be 
written as 



^Tr In V* (V) = ^Tr - W v - m 2 ^) - '-Tr (m 2 ^ - V„V*) 



(13) 



We can see in ({TBI) that the first term is the effective action of a minimal second order oper- 



ator K£ (V) with potential —R%. The second term can be transformed as Tr 



Tr 



■iVD n - 1 V l 



= Tr 




n 

and 






-Tr (' 
2 v 


-r?X - V„V") = ^Tr (m 2 v - □) . 


(14) 



Then, the effective action for the massive vector field is equal to the effective action of the 
minimal second order operator K£ (V) minus the effective action of a minimal operator 
&v (V) corresponding to a massive scalar field minimally coupled to gravity. 

The problem with the Dirac nonminimal operator Df is solved introducing a new spinor 
variable ip connected with <p by the relation <p = 7°"V CT ?/> — m/0 so that (jSJ) take the form 
T'V V m V„V> - mji/; = 0. Making use of the identity 7^7" V^V*, = I (□ - \R) equation (J8J 
becomes of the form 



f\rnin 



□ 



m, 



if) = 



(15) 



\RI. 



where the potential matrix can be easily identified as Q - 

Now using the Schwinger-DeWitt representation for the Green functions of the minimal 
operators, we can obtain for the renormalized one loop effective action of the quantum mas- 



sive vector and spinor fields the expression T^ ren 
effective Lagrangian reads: 



-g 2, re n where the renormalized 



2(4*) 2 k(k-l)(k -2) 



Tra^^XjX) — Tr a^ u; (a;, x) ^ Tr a K k 2j (x, x) 



] = a^(x,z'), [a^ 1 ] = a^'(x,x') and [a^] 



,(°) 



(0), 



2(fc-2) 



m 



2(fc-2) 
/ 



, (16) 



.(§)/ 



whose coincidence limit 



appears under the supertrace operation in ([16]) are the HMDS coefficients for the minimal 
operators K, S and IW in respectively. As usual, the first three coefficients of the DeWitt- 
Schwinger expansion, a , a 1; anda2, contribute to the divergent part of the action and can 
be absorbed in the classical gravitational action by renormalization of the bare gravitational 
and cosmological constants. 

Restricting ourselves here to the terms proportional to m ~ 2 , using integration by parts 
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26(|, we obtain for the 



and the elementary properties of the Riemann tensor 
renormalized effective lagrangian in the case of the massive vector field considered in this 
work 
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(17) 



The interested reader can find the general result for the spinor field case, for example, in 



reference 



171 ] . As we can see, this final expression of the one loop effective for the massive 



vector field only differ from that of the massive scalar and spinor fields in the numerical 
coefficients in front of the purely geometric terms. For (T^ u ) ren we obtain a very cumbersome 
expression that, as in the case of ( fl7j) . is different from that obtained for scalar and spinor 
fields only in the numerical coefficients that appears in front of the purely geometrical terms. 
For this reason we not put this very long expression for the stress tensor here and refers the 



readers to our previous papers 
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For the present work we deal with the Reissner-Nordstrom spacetime. The obtained 
results for the components of the stress-tensor are very simple and can be found in reference 
18]. 



For a quantum scalar field (f>(x) with mass m interacting with gravity with non minimal 



coupling constant £ we find 

{T.K = CS+(t-^DS , (18) 

where 

C\ = -T S (1248Q 6 - 810r 4 Q 2 + 855M 2 r 4 + 202r 2 Q 4 - 1878M 3 r 3 + 1152Mr 3 Q 2 
+ 2307M 2 r 2 Q 2 - 3084MQ 4 r), 

L>* = S(-792M 3 r 3 + 360M 2 r 4 + 2604M 2 Q 2 r 2 - 1008MQV - 2712MQ 4 
+ 819Q 6 + 728Q 4 r 2 ), 
C; = T S (444Q 6 - 1488MQV + 162Q 2 r 4 + 842Q 4 r 2 - 1932MQ 4 r + 315M 2 r 4 
+ 2127M 2 Q 2 r 2 - 462M 3 r 3 ), 
D r r = S(-792M 3 r 3 + 360M 2 r 4 + 2604M 2 Q 2 r 2 - 1008MQV - 2712MQ 4 r 

+ 819Q 6 + 728Q 4 r 2 ), 
C e e = -T S (3044QV - 2202M 3 r 3 - 10356MQ 4 r + 3066Q 6 - 4884MQV 

+ 9909M 2 Q 2 + 945M 2 r 4 + 486Q 2 r 4 ), 
D% = 2(3276M 2 Q 2 r 2 - 1176MQ 2 r 3 - 3408MQ 4 r + 1053Q 6 - 1008M 3 r 3 
+ 432M 2 r 4 + 832Q 4 r 2 ), 
For the massive vector field we obtain: 

(T v )l = -T 1) (-31057Q 6 - 12150r 4 Q 2 - 1665M 2 r 4 - 41854r 2 Q 4 - 93537M 2 r 2 Q 2 

+3666M 3 r 3 + 69024Mr 3 Q 2 + 10751MQ 4 r) , 
(T v ) r r = T 1) (-10448MQ 2 r 3 + 2430QV + 6442QV - 693M 2 r 4 + 12907M 2 Q 2 r 2 

+5365Q 6 - 16996MQ 4 r + 1050M 3 r 3 ) , 
(T v ) e e = -T^(20908Q 4 r 2 + 4854M 3 r 3 - 44068MQ 4 r - 31708MQ 2 r 3 + 7290Q 2 r 4 
+30881M 2 Q 2 r 2 - 2079M 2 r 4 + 13979Q 6 ) , 
and for the spinor field the resulting components are given by 

(T f )l = -T / (-21496 rMQ 4 + 4917 Q 6 + 10544 r 2 Q 4 - 22464 Mr 3 Q 2 + 21832 M 2 r 2 Q 2 

-1080 M 2 r 4 + 2384 M 3 r 3 + 5400 r 4 Q 2 ) , 
(Tf) r r = T7(-6336 Mr 3 Q 2 + 8440 M 2 r 2 Q 2 + 2253 Q 6 + 3560 r 2 Q 4 - 8680 rMQ 4 
+504 M 2 r 4 - 784 M 3 r 3 + 1080 r 4 Q 2 ) , 
(T f } g = -17(12080 r 2 Q 4 - 33984 rMe 4 + 9933 Q 6 + 30808 M 2 r 2 Q 2 - 20016 Mr 3 Q 2 
+ 1512 M 2 r 4 - 3536 M 3 r 3 + 3240 r 4 Q 2 ) . 
In the above expressions, we have T s = (302407r 2 m 2 r 12 )™ 1 ,T V = (100807r 2 m 2 r 12 )~ 1 ,T/ = 

8 



(403207r 2 m^r 12 ) and 5 = (720tt 2 m 2 r 1 2) . Q and M denotes the charge and bare mass 
of the black hole. 

The above quantum stress energy tensors are regular at the event horizon, as is to be 
expected due to the local nature of the Schwinger-DeWitt approximation and the regular 
nature of the horizon. Also from the general form of the geometric terms conforming the 
general expresion for the constructed stress tensor, we see that it is covariantly conserved, 
thus indicating that it is a good candidate for the expected exact one in our large mass 
approximation. 

In the specific case of a Reissner-Nordstrom spacetimes, Anderson et.al showed, using 
detailed numerical results for the scalar field case, that for m s M > 2 the deviation of the 
approximate stress energy tensor from the exact one lies within a few percent [ijj]. The 
more general condition for the validity of the Schwinger-DeWitt approximation in the case 
of a spin-j field can be written as rrijM > 1, where rnj and M are respectively the field and 
black hole masses. In the following we carefully take into account the above condition in 
numerical calculations. 



III. SEMICLASSICAL SOLUTION 



In our previous paper we have shown how to find the general solution to the backreaction 



problem in spacetimes with spherical symmetry [lOj. In the following we solve the general 
semiclassical Einstein equations assuming that there are an electromagnetic field as a clas- 
sical source, and multiple quantized massive free fields as a perturbative quantum source, 
so the solution to the backreaction problem gives a quantum corrected Reissner-Nordstrom 
black hole. 

In this limit in which we deal with free fields on a background spacetime, the treatment 
of the backreaction due to a collection of fields with different spins is easy. This is due to the 
fact that the quantum stress tensors in this limit only depends quadratically on the fields 
and are flavour diagonal. For example, for a set of N s real scalars, upon renormalization we 
have 



v 



< r ;>L = E(( T ;) fc L = ^<( T ;)>L > ( 19 ) 
fc=i 

where (TJZ) is the classical stress tensor for a single scalar. The last equality above follows 



9 



from the fact that the renormalization procedure is independent of the species label k. In a 
similar manner we can arrive to the same conclusion for the renormalization of spinor and 
vector field stress tensors. 

The above statements permit us to obtain a good approximation to the multiple field 
backreaction using as the source term in the semiclassical Einstein equations an appropriate 
weighted combination ( with weights N s , N v and Nf ) of the single-species renormalized 
stress energy tensors. 

In the case of our interest the general form for the line element that solves the backreaction 
problem, considering only terms that are linear in the perturbation parameter e = 1/M 2 , is 
given by 

ds 2 = -A{r)dt 2 + B{r)dr 2 + r 2 (d6 2 + sin 2 6 d</) 2 ) , (20) 



with 



and 



where 



B(r) r 



^4 + 7^f('(^ ■ <*) 

3 J °° 



A ( r ) = ^II ex P^( r )} > ( 22 ) 



A,-(r) = SvriV, f (B (() (( - ( 7?) ) d( . (23) 

Joe 

where the subindex j denotes the different single spin species considered (scalar, vector and 
spinor field). Inserting the corresponding expressions for the temporal component of the 
quantum stress tensor for the different fields considered in this work, we obtain 

1 2M Q 2 e JV-- , . 

-sn = 1 + -r + -E^> r , 24 

B(r) r r l n ^-^ mi 

j j 

where, for the scalar case 

F s (r) = E(r)+ZH(r) , (25) 

with 

613M 3 Q 4 2327M 2 Q 6 3M 2 Q 2 5M 4 1237M 5 883M 2 Q 4 

^ ~ 840r 9 + 1134r 10 70r 6 + 28^ ~ 3780r 7 + 4410r 8 

82M 3 Q 2 1369M 4 Q 2 

315r 7 + 1764r 8 ' 
_ 28M 3 Q 2 113M 3 Q 4 91M 2 Q 6 52M 2 Q 4 4M 4 22M 5 62M 4 Q 2 

^ ~~ 15r 7 + 30r9 90^0 ior 8 5r^ + 15r 7 15r 8 ' 
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e /184M 3 Q 2 29M 4 229M 2 Q 4 M 2 Q 2 

s ^ ~~ T^ 2 V 441r 7 140r 6 840r 8 + 35r 6 

e£ /14M 4 13M 2 g 4 32M 3 Q 2 \ 

7rm 2 \ 15r 6 10r 8 15r 7 / ' 

and for vector and spinor fields 

_ 26879 M 3 g 4 2876 M 3 Q 2 611 M 5 37 M 4 20927 M 2 Q 4 

~ 2520 r 9 + "315" r 7 + 1260 T 7 " ~ 140 "r^ ~ 4410 r 8 
10393 M 4 Q 2 27 M 2 Q 2 31057 M 2 Q 6 

~ ""980 r 8 14 ""^ 11340 r 10 ' 
2687 M 3 Q 4 1639 M 2 Q 6 149 M 5 

/ ^ ~ 5040 15120 ~r To 1890 ~ 

3 M 2 Q 2 3 M 4 26 M 3 Q 2 2729 M 4 Q 2 659 M 2 Q 4 



14 r 6 70 r 6 35 r 7 4410 r 8 2205 r 8 

e / 131 M 4 9 M 2 Q 2 _ 9784 M 3 Q 2 2141 M 2 Q 4 \ 

v ^ ~ nm v 2 V420 7^ + 7 "7^ 2205 r 7 + W r 8 J 
e / 37 M 2 g 4 11 M 4 1 M 2 Q 2 52 M 3 Q 



(26) 



7rm f 2 V560 r 8 210 r 6 + 7 r 6 245 r 7 
The horizon for the quantum corrected solution will be, up to first order in the perturbation 

parameter, at position r + given by 

(27) 



where 



(M - Q 2 /r H ) 



rh 



^ = - 7i7Z^3 / C*<I?(0>,««C ■ (») 



and r# is the position of the event horizon of the bare classical solution. Upon substitution 
of the required quantities in the above expression we find 

£ Y ■ 

Aj = nm 2 (M-Q 2 /r H ) ' (29) 

with 



v s = e + f n, (30) 

with 

_ 613M 3 Q 4 2327M 2 Q 6 3M 2 Q 2 5M 4 1237M 5 883M 2 Q 4 
1680r| 22680r^ + U0r 5 H ~ 5Qr% + 75Q0r% ~ 8820r^ 
41M 3 Q 2 1369M 4 Q 2 
+ ~315r| 3528^ ' 

14M 3 Q 2 113M 3 Q 4 91M 2 Q 6 26M 2 Q 4 2M 4 11M 5 31M 4 Q 2 
I5r% 60r 8 + 180r 9 + 45rL + ~ 15r 6 w + 15r 7 
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and for the vector and fermion components we find 

_ 27 M 2 Q 2 37 M 4 26879 M 3 Q 4 1438 M 3 Q 2 31057 M 2 Q 6 

v ~ 28 r^ 5 + 280 7^ ~ 5040 r^ 8 ~ "315" r H e + 22680 r# 9 

611 M 5 20927 M 2 Q 4 10393 M 4 Q 2 



2687 M 3 Q 4 149 M 5 3 M 4 3 M 2 Q 2 2729 M 4 Q 2 
r f = -tt^t: T~ + tt^tt — = - ttt: — I + ™ a + 



(31) 



10080 r H 7 3780 r H 5 140 r^ 4 28 r H 4 8820 r H e 
13 M 3 Q 2 659 M 2 Q 4 1639 M 2 Q 6 
35 r H 5 + 4410 r^ 6 + 30240 r H 8 



IV. SCALAR PERTURBATIONS AND QUASINORMAL MODES 



(32) 



In the following we consider the evolution of a test massless scalar field $(x M ) with 

= (t,r,0,(j>), in the quantum corrected gravitational background studied above. The 
dynamics of for this test field is governed by the Klein-Gordon equation 

^rZ; 1=0 , (33) 



with is the metric tensor of semiclassical solution and g its determinant. Upon separation 
of the angular and radial part in the above equation and the introduction of the radial tortoise 
coordinate 

£jZ L -[ U *-V]Z L = , (34) 

where Zl(t) denotes the radial component of the wave function, oj is the quasinormal fre- 
quency and V is the effective potential. The potential V is a function of the metric compo- 
nents g^ u and the multipolar number L, and for the test massless scalar field considered in 
this work, is given by the general expression 

V[r(r.)]=A(r)y^ + ^^ , (35) 

where the prime refers to the derivative with respect to the radial coordinate r. For semi- 
classical black holes we have in general the following expression for the effective potential 

V{r) = V c {r) + Uj{r) + O (e 2 ) , (36) 

where V c (r) is the scalar effective potential of the bare black hole solution and U(r) is a 



complicated function of the contributions -^Ujir) due to the vacuum polarization effect 
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related the multiple field backreaction, that we will not write here. In the case of a classical 
Reissner-Nordstrom black hole V c (r) is given by 



r 



2 



2 Mr + Q 2 ) (-2 Q 2 + (3 r 2 + 2 Mr) 



„c (r)= v " 7V r6 " ^ , (37) 

where (5 = L(L + 1). For the semiclassical black hole solution considered in this paper, 
where the vacuum polarization effects comes from the quantization of massive scalar, vector 
and spinor fields in the large mass limit, the particular expression for Uj(r) results 

U a (r) = W 1 (r)+ZW 2 (r), (38) 

where 



WAr) 



+ 



1751 M 2 Q 4 9 M 4 1021 M 5 1816 M 6 6 M 2 Q 2 

4410 ~Y° 20 ~Y + "Yo" ~Y ~ "945" Y^ + 35 r 8 
674641 M 3 Q 6 17 M 3 Q 2 13271 M 4 Q 4 625 M 2 Q 8 



158760 r 13 105 r 9 1764 r 12 756 r 14 

23353 M 2 Q 6 8559 M 3 Q 4 962 M 4 Q 2 16687 M 5 Q 2 

~ 15876 r 12 + I960 Y 11 245 r 10 + 2940 r 11 

/ 1529 M 2 Q 6 1 M 4 1 M 2 Q 2 
+ L [L + lj I — 



22680 r 12 35 r 8 70 r 8 
47 M 5 773 M 2 Q 4 44 M 3 Q 2 821 M 3 Q 4 1171 M 4 Q< 



and 



540 r 9 17640 r 10 441 r 9 3528 r 11 4410 r 10 

4 M 3 Q 2 16 M 4 Q 2 29 M 3 Q 4 2 M 5 13 M 2 Q 4 



15 r 9 15 r 10 30 r 11 5 r 9 90 r 10 
13 M 2 Q 6 2 M 4 \ 26 M 2 Q 6 162 M 5 Q 2 2101 M 3 Q 6 



45 r 12 15 r 8 y 3 r 12 5 r 11 90 r 13 
128 M 6 M 4 13 M 2 Q 8 182 M 2 Q 4 289 M 3 Q 4 28 M 3 Q 2 
15 r 10 r 8 3 r 14 45 r 10 10 r 11 5 r 9 



42 M 5 416 M 4 Q 2 130 M 4 Q 



4 



+ -rz + 



5 r 9 15 r 10 3 r 12 
for the scalar case and 

107577 M 5 Q 2 306442 M 4 Q 2 34907 M 4 Q 4 13 M 4 54 M 2 Q 2 

v ^ ~ ~ ~980 + 2205 r 10 + 196 r 12 + 20 "Y" + Y r 8 

491 M 5 872 M 6 1205 M 3 Q 2 303071 M 3 Q 4 15415961 M 3 Q 6 



180 r 9 315 r 10 21 r 9 1960 r 11 158760 r 13 

681461 M 2 Q 6 29027 M 2 Q 4 66421 M 2 Q 8 / 5M 5 

+ 15876 Y 42 "" + "882" "Y^" + "378Y ~r Tr ~ + ( + l > \ 36 Y~ 

4678 M 3 Q 2 6653 M 3 Q 4 16067 M 2 Q 4 4307 M 2 Q 6 1 M 4 

+ 2205 r 9 + 5880 Y 11 17640 ~Y° 22680 r 12 + 21 ~Y 
6257 M 4 Q 2 9 M 2 Q 2 ' 

~ 4410 ~Y° 14 r 8 



(39) 
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and 



Uf(r) = - 



2279 M 2 Q 8 11101 M 5 Q 2 24764 M 4 Q 2 5092 M 4 Q 4 191 M 3 Q' 



2520 r 14 

'4 



1470 r 11 2205 r 10 441 r 12 35 

r3/o6 iiq *r5 a i\/r2/^,2 



1 M 4 119141 M 3 Q e 113 M 5 6 M' 2 Q' 2 80 M e 8775 M 3 Q 4 
10 T 8 " ~ 21168 r 13 + 270 7^ + 7 r 8 



23797 M 2 Q 4 3569 M 2 Q 6 
+ "8820" ~7^ + 1323 "7^ + L ( L + 1 } 
1 M 2 Q 2 3173 M 2 Q 4 544 M 4 Q 2 



14 r s 



35280 r 



10 



2205 r 



10 



r 8 189 r 10 784 r 11 
7 M 5 1 M 4 12 M 3 g 2 

270 7^ ~ 105 78" + 49 r 9 

8 M 2 Q 6 6659 M 3 Q 4 



189 r 12 



35280 r 



ii 



(40) 



for the vector and spinor cases. 

In Figure ([[]) is presented the effective potential V taking into account, as an example, 
the backreaction of multiple fields with weights N s , N v and Nf all equal to 10. 



Vef f [r; 

-a 

3-10 



2.5-10" 
2 -10" 

1.5-10" 
1 - 10" 
5 - 10" 



L = 




200 400 600 800 1000 1200 



r 



Figure 1. Effective potential of L=0 scalar modes for Semiclassical black hole with M = 100 (top 
curve) and M = 200(bottom curve). Q/M = 0.95, N s = N v = Nf = 10 and the mass parameter 
of all quantum fields are chosen to be m=l/10. 



As it is observed, the figure shows a definite positive potential barrier, i.e, a well behaved 
function that goes to zero at spatial infinity and gets a maximum value near the event 
horizon. The quasinormal modes are solutions of the wave equation ( 1341) with the specific 
boundary conditions requiring pure out-going waves at spatial infinity and pure in-coming 
waves on the event horizon. The quasinormal frequencies are in general complex numbers, 
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whose real part determines the real oscillation frequency and the imaginary part determines 
the damping rate of the quasinormal mode. In order to evaluate the quasinormal modes 
for the situation considered in this report field, we use the well known WKB technique at 
sixth order, that can give accurate values of the lowest ( that is longer lived ) quasinormal 
frequencies, and was used in several papers for the determination of quasinormal frequencies 
in a variety of systems j^ . The first order WKB technique was applied to finding quasi- 



normal modes for the first time by Shutz and Will 



Latter this approach was extended 



to the third order beyond the eikonal approximation by Iyer and Will 29[ and to the sixth 



order by Konoplya 30]. 



The sixth order WKB expansion gives a relative error which is about two orders less than 
the third WKB order, and allows us to determine the quasinormal frequencies through the 
formula 

i ^^.-±U j =n+\ , , = 0,1,2,... (41) 

where Vq is the value of the potential at its maximum as a function of the tortoise coordinate, 
and Vq represents the second derivative of the potential with respect to the tortoise coor- 
dinate at its peak. The correction terms IL depend on the value of the effective potential 



and its derivatives ( up to the 2i-th order) in the maximum, see 3JJ and references therein. 

The results of the numerical evaluation of the first two fundamental quasinormal frequen- 
cies in the case considered in this work is showed in table (JT]). 

From figures ([2]) and ([3]) we see that the backreaction of the quantized multiple fields upon 
the classical charged black hole gives rise to an increasing of the real oscillation frequencies 
and to a decreasing of the damping rate, for physically interesting values of the black hole 
mass. Then, as a consequence of the vacuum polarization effect due to the multiple massive 
fields, we have an effective increasing of the quality factor, proportional to the ratio r^fe^i ■ 
As expected, the differences in the quasinormal frequencies when the black hole mass in- 



creases tend to become small. It is interesting to note that the above results show significant 
differences with those obtained previously considering only the backreaction of a quantized 
massive scalar field upon the classical charged black hole solution. In that case the quality 
factor is small for semiclassical black hole with respect to the classical solution, and the shift 
in the quasinormal frequencies is less pronounced. As the numerical calculations show, this 
is not the case if the separate backreaction due to massive vector and spinor fields are con- 
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Semiclassical solution 


Classical solution 








M = 


100 






L 


n 


Re(w) 


—Im{w) 


L 


n 


Re(w) 


—Im{w) 








10.2421 


6.64263 








5.52131 


29.4487 


1 





28.7831 


6.80833 


1 





27.1215 


7.5690 


1 


1 


26.5433 


21.0202 


1 


1 


23.2393 


33.6626 








M = 


120 






L 


n 


Re(w) 


—Im{w) 


L 


n 


Re(w) 


—Im{w) 








8.53509 


5.53553 








4.6011 


24.5404 


1 





23.9859 


5.6736 


1 





22.2679 


6.30758 


1 


1 


22.1194 


17.0889 


1 


1 


19.3661 


28.0521 








M = 


150 






L 


n 


Re(w) 


—Im{w) 


L 


n 


Re(w) 


—Im{w) 








6.82809 


4.42843 








3.68089 


19.6324 


1 





19.1887 


4.5389 


1 





19.0144 


5.04608 


1 


1 


17.6956 


14.0135 


1 


1 


15.4929 


22.4418 



Table I. Rescaled scalar quasinormal frequencies w = 10 w for the classical and semiclassical 
Reissner-Nordstrom black hole, with Q/M = 0.95, N s = N v = Nf = 10 and m = 1/10. 



sidered. In each one of this cases, we obtain an increasing of the quality factor more relevant 
that the decreasing showed by the scalar field case. As a result, the combination of all types 
of fields give a higher quality factor with respect to the bare black hole. Thus, we arrive to 
the conclusion that the semiclassical Reissner-Nordstrom black holes are better oscillators 
than its classical partners. Then, in this case, there is a qualitative correspondence with the 



results obtained by Konoplya in reference ll| for the BTZ black hole dressed by a quantum 



conformal massless scalar field, where he studied the backreaction due to Hawking radiation 
upon the classical spacetime. 

We also studied the dependence of the quasinormal frequencies for a fixed black hole bare 
mass and different values of the the quantum field mass rrij, obtaining similar results with 
respect to the massive scalar field case: a little dependance of the quasinormal frequencies 
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Figure 2. Dependance of Re{oj) on M for classical and semiclassical black holes. The parameters 
are chosen to be Q/M = 0.95, m = 1/10, L = 0, N s = N v = N f = 10 and n = 0. 



on this parameters. As the quantum field mass increases, we found a very small increment 
in the real part of the frequencies for semiclassical black holes, and a very small decrease in 
the imaginary part. The same occur if we consider the variation of the coupling constant 
between the massive scalar field component and the gravitational background: the quasinor- 
mal frequencies are insensitive to the variation of this parameter. Therefore, the shift in the 
quasinormal frequencies with respect to the classical bare black hole appears to be the same 
for the given range of the quantum field masses. With respect to the multiplicity number Nj 
of a given quantum field we find that, as expected, the shift shows some increment as this 
parameters increases, an effect that is more pronounced for very large Nj. It is important 
to take care with the fact that, very large values of this numbers can imply that the total 
quantum stress tensor becomes a large quantity, such that the perturbative treatment of the 
backreaction problem used here becomes inadequate. By fortune, this is not the case for 
physically interesting values of Nj. 
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Figure 3. Dependance of Im(uj) on M for classical and semiclassical black holes. The parameters 
are chosen to be Q/M = 0.95, m = 1/10, L = 0, N s = N v = N f = 10 and n = 0. 



V. CONCLUDING REMARKS 



We have studied the influence of the backreaction due to vacuum polarization of multiple 
species of large mass quantum massive fields, belonging to the standard model, upon the 
structure of scalar quasinormal frequencies for semiclassical charged black holes. The effect 
observed is a shift in the quasinormal frequencies the semiclassical solution, such that quan- 
tum corrected black holes becomes better oscillators that classical ones. It is important to 
verify that the above effects are also true for the quasinormal ringing phase in the evolution 
of spinor and electromagnetic test fields. We are currently investigating such problems and 
the results will be presented in future reports. 
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